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Abstract 

Steady linear gravity waves of small amplitude travelling on a current of constant 
vorticity are found. For negative vorticity we show the appearance of internal waves 
and vortices, wherein the particle trajectories are not any more closed ellipses. For 
positive vorticity the situation resembles that of Stokes waves, but for large vorticity 
the trajectories are affected. 
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1 Introduction 

The subject of this paper are periodic gravity water waves travelling with 
constant shape and speed. Such wave-trains are an everyday observation and, 
typically, one gets the impression that the water is moving along with the 
wave. In general, this is not so. Rather, the individual fluid particles display a 
motion quite different from that of the wave itself. While for irrotational waves, 
recent studies have enlightened the situation, we investigate the situation of 
waves propagating on a rotational current, so that there is a non- vanishing 
curl within the velocity field. 
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For irrotational waves, there is a classical first approximation showing that the 
fluid particles move in ellipses, back and forth as the wave propagates above 
them. This can be found in classical [1,2,3] as well as modern [4,5] text books, 
and it is consistent with the only known explicit solutions for gravity water 
waves: the Gerstner wave [6,7] for deep water, and the edge wave solution for a 
flat beach [8] , both with a depth- varying vorticity. A formal physical argument 
involving a balance between opposing forces was used in [9] to get a similar 
result without the use of irrotationality. There arc also experimental evidence 
supporting this picture. Those include photographs [4,10,3] and movie films 
[11]. 

However, as anyone having used bottle post would guess, there are other find- 
ings. Even in [9], where it is asserted that the orbits are elliptic, and where the 
photographs and movie films are referenced, the author notes that "/ am not 
aware of any measurements that show that the particle orbits of shallow water 
waves are indeed ellipses." In fact it was observed already in the 19th century 
that there seems to be a forward mass drift [12], so that the average motion 
of an average fluid particle is along with the wave. This phenomenon can be 
seen by making a second approximation of the governing equations, and it is 
known as Stokes drift (see also [13,14]). In [15,16] it was deduced that for steep 
waves the orbits deviate from simple ellipses. There is also mathematical evi- 
dence uniformly showing that a more thorough study of the equations yields 
non-closed orbits with a slight forward drift. Those include investigations of 
the precise orbits of the linearized system [17,18,19], and two recent papers 
on exact Stokes waves [20,21] (steady irrotational and periodic gravity waves 
which are symmetric and monotone between trough and crest). The relation 
between such results and experimental data is discussed in [20], where it is 
argued that the ellipses - at least near the bottom - are approximations of 
the exact trajectories. 

While many situations are adequately modelled by irrotational flows - e.g. 
waves propagating into still water - there are situations when such a math- 
ematical model is insufficient. Tidal flow is a well-known example when con- 
stant vorticity is an appropriate model [22], a fact conflrmed by experimental 
studies [23]. This is one reason why, recently, the interest for exact water 
waves with vorticity has increased. At this point existence [24], variational 
characterization [25], uniqueness [26,27], symmetry [28], and a unique contin- 
uation principle [29] for flnite depth steady gravity waves with vorticity are 
established. There is also a theory for deep-water waves [30,31], as well as for 
capillary and capillary-gravity waves [32,33]. However, due to the intricacy of 
the problem, studies of the governing equations for water waves are extremely 
difficult. In-depth analyses are very rare. To gain insight into qualitative fea- 
tures of flows with vorticity Ko and Strauss recently performed a numerical 
study [34] extending earlier work by DaSilva and Peregrine. We will pursue a 
different approach. Notice that the intuitive notion of vorticity is captured in 
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what happens when one pulls the tap out of a bath tub. It should therefore 
come as no surprize that the particle paths of waves travelling upon a rota- 
tional current deviate from those in the case of waves without vorticity. That 
is the main result of this paper. More precisely, we make a first attempt at 
understanding the particle trajectories by deducing a linear system for con- 
stant vorticity. Here, linearity means that the waves arc small perturbations 
of shear flows, hence of small amplitude. The system obtained is solvable in 
the sense of closed expressions, and thus it is possible to make a phase portrait 
study of the steady wave. 

It is found that for positive vorticity, the steady wave resembles that of the 
irrotational situation [17,18], though the physical particle paths behave differ- 
ently if the size of the vorticity is large enough. For negative vorticity, however, 
we show the existence of a steady periodic surface wave containing an internal 
wave as well as a vortex, or so called cats-eye (cf. [35] and [36, Ex. 2.4]). For 
unit depth this situation occurs if the absolute size of the negative vorticity 
exceeds the wave speed, while in the opposite situation both the steady wave 
and the physical particle trajectories resemble the irrotational case. When the 
size of the negative vorticity exceeds the wave speed the particle trajectories 
of the internal wave behave in the same manner as in the irrotational case - 
nearly closed ellipses with a forward drift - but within the vortex and the sur- 
face wave the particles are moving mainly forward. This indicates that such a 
wave may be unordinary or unstable, since measurements show that for waves 
not near breaking or spilling the speed of an individual particle is generally 
considerably less than that of the wave itself [2] . Such a situation is excluded 
in [20,21], and our result is therefore not in contrast to those investigations. 

An interesting feature of the phase portrait for negative vorticity is that it 
captures the almost ideal picture of what vorticity is. It furthermore indicates 
that in the case of large negative vorticity the governing equations allow for 
travelling waves very different from the classical Stokes waves (see [37] for 
a good reference of that subject). Finding those waves with analytic tools 
could prove difficult; so far the existence results [24,31,38] for steady waves 
with vorticity rely on the assumption that no particle moves as fast as the 
wave itself. This study suggests that the presence of vorticity - even when it 
is constant - changes the particle trajectories in a qualitative way, that this 
change depends on the size of the vorticity, and that it applies less to particles 
near the bottom. 

The disposition is as follows. Section 2 gives the mathematical background for 

the water wave problem, while in Section 3 we deduce the linearization and its 
solution. The main findings are presented in Section 4, and the implications 
for the particle trajectories in Section 5. In Section 6 we give a brief summary 
and discussion of our results. 
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2 Preliminciries 



The waves that one typically sees propagating on the surface of the sea are 
locally approximately periodic and two-dimensional (that is, the motion is 
identical in any direction parallel to the crest line). Therefore - for a descrip- 
tion of these waves propagating over a flat bed - it suffices to consider a cross 
section of the flow that is perpendicular to the crest line. Choose Cartesian 
coordinates {x, y) with the ?/-axis pointing vertically upwards and the x-axis 
being the direction of wave propagation, while the origin lies on the flat bed 
below the crest. Let (u(t,x,y), v(t,x,y)) be the velocity fleld of the flow, let 
h > he the depth below the mean water level y = h, and let y = h + rjit, x) 
be the water's free surface. We assume that gravity is the restoring force once 
a disturbance was created, neglecting the effects of surface tension. Homo- 
geneity (constant density) is a physically reasonable assumption for gravity 
waves [2], and it implies the equation of mass conservation 

Ux + Vy — (2.1a) 

throughout the fluid. Appropriate for gravity waves is the assumption of in- 
viscid flow [2] , so that the equation of motion is Euler's equation 

ut + uu^ + vuy = -Px, ^2 i^-^ 

Vt + uv^ + vvy = -Py - g, 

where P{t, x, y) denotes the pressure and g is the gravitational constant of 
acceleration. The free surface decouples the motion of the water from that of 
the air so that, ignoring surface tension, the dynamic boundary condition 

P = Po on y = h + 7j{t,x), (2.1c) 

must hold, where Pq is the constant atmospheric pressure [5] . Moreover, since 
the same particles always form the free surface, we have the kinematic bound- 
ary condition 

v^r)t + ur)^ on y = r]{t, x). (2- Id) 

The fact that water cannot penetrate the rigid bed at y = yields the kine- 
matic boundary condition 

v = on y = 0. (2.1e) 
The vorticity, cu, of the flow is captured by the curl, 

Vx — Uy — LO. (2-lf) 



We now introduce a non-dimensionalization of the variables. As above, h is the 
average height above the bottom, and we let a denote the typical amplitude, 
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and A the typical wavelength. It is reasonable - and fruitful - to take \/gh. as the 
scale of the horizontal velocity. That is the approximate speed of irrotational 
long waves [5]. We shall use c to denote the wave speed, and we let 



gh 

be the starting point of the non-dimensionalization. We then make the trans- 
formations 

X y \fght u \v rj 

Ah A ^ygl^ n^gn a 

Having made these transformations, define furthermore a new pressure func- 
tion p — p{t, X, y) by the equality 

P = Po + 9h{l -y)+ ghp. 

Here Pq is the constant atmospheric pressure, and gh{l — y) is the hydrostatic 
pressure distribution, describing the pressure change within a stationary fluid. 
The new variable p thus measures the pressure perturbation induced by a 
passing wave. It turns out that the natural scale for the vorticity is ^Jh/g and 
we thus map 

(h 

a; I— > 4/ — a;. 
V 9 

The water wave problem (2.1) then transforms into the equations 

u,^ + Vy = 0, (2.2a) 
Ut + uu^ + vuy = -p^, (2.2b) 



A^ 

vt + uvx vvy = ~j^Py^ (2.2c) 
h 



^^v^-Uy=uj, (2. 2d) 



valid in the fluid domain 0<y< H-^ry, and 

V ^ ^iVt + UTj^), (2.2e) 
P^lv, (2.2f) 

valid at the surface y = 1 + |?7, in conjunction with the boundary condition 
(2.1e) on the flat bed y = 0. Here appear naturally the parameters 

a . h 

called the amplitude parameter, and the shallowness parameter, respectively. 
Since the shallowness parameter is a measure of the length of the wave com- 
pared to the depth, small 5 models long waves or, equivalently, shallow water 
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waves. The amplitude parameter measures the relative size of the wave, so 
small e is customarily used to model a small disturbance of the underlying 
flow. We now set out to study steady (travelling) waves, and will therefore 
assume that the equations (2.1) have a space-time dependence of the form 
X — ct va. the original variables, corresponding to X{x — ct) in the equations 
(2.2). The change of variables 



yields the problem 



{x,y) ^ {x- ct,y) 

Ux + Vy = 0, (2.3a) 

{u - c)ux + vuy = -px, (2.3b) 

{u- c)vx + vvy = (2.3c) 

S'^Vx — Uy = LU, (2.3d) 



valid in the fluid domain < y < 1 + er), 



V = e{u- c)r]^, (2.3e) 
P = ev, (2.3f) 



valid at the surface y — 1 + erj, and 



along the flat bed y — 0. 



(2.3g) 



3 The linearization 



To enable the study of explicit solutions, we shall linearize around a laminar 
- though rotational - flow. Such shear flows are characterized by the flat 
surface, y = 1, corresponding to ?7 = 0, so insertion of this into (2.3) yields 
the one-parameter family of solutions, 

ry 

U{y) = U{y;s) = s - uj{y)dy, 

Jo 

with T] = 0, p = 0, V = 0. We now write a general solution as a perturbation 
of such a solution U, i.e. 

u — U + eu, V — ev, p — ep. (3.1) 

We know from the exact theory of water waves that such solutions exist at 
the points where the non-trivial solutions bifurcate from the curve of trivial 
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flows [24]. Remember that small e corresponds to waves whose amplitude is 
small in comparison with the depth. Since the surface is described by 1 + er], 
Tj should thus be of unit size. Dropping the tildes, we obtain 

Ux-\-Vy — 0, (3.2a) 
{U — c)ux + vUy + e{vuy + uux) — —Px, (3.2b) 

{U - c)vx + e{vvy + uvx) ^ (3.2c) 
valid in the fluid domain < y < 1 + erj, 

V = {U -c + eu)rix, (3.2d) 
P = V, (3.2e) 

valid at the surface y — 1 + er), and 

v = (3.2f) 

on the flat bed y = 0. The corresponding linearized problem is valid in the 
sense that its solution satisfies the exact equations except for an error whose 
size can be expressed as a square of the size of the linear solution. The lin- 
earization is attained by formally letting £ — > 0, and it is given by 

Ux + Vy — 0, (3.3a) 
{U - c)Ux + vUy = -px, (3.3b) 

{U-c)vx^-^, (3.3c) 



valid for < y < 1, and 



v^{U- c)r]x, (3.3d) 
p^rj, (3.3e) 

valid for y = 1. In order to explicitly solve this problem we restrict ourselves 
to the simplest possible class of vorticities, i.e. when Lj{y) = a; G M is constant. 
It then follows that 

U{y;s) = -(jy + s. 

Looking for separable solutions we make the ansatz r}{x) — cos {27rx) (note 
that the original wavelength A and the original amplitude a have both been 
non-dimensionalized to unit length). The solution of (3.3) is then given by 

'u{x,y) = 2SnC cos {2'7ix) cosh {27i6y) , 
v{x, y) = 2ttC sin {2-kx) sinh {2Ti5y) , 

\p{x^ y) — C cos {2txx) {2'k8{c — s + ujy) cosh {27rSy) — uj sinh {27rSy) ^ , 

(3.4) 



7 



where 

^ _ c - S + UJ 
sinh(27r(5) ' 
and c, 5, /i, s, uj must satisfy the relation 

{c-s + uj) {2ti5{c -s + uj) coth(27r5) - uj) = I (3.5) 

This indicates that the properties of the wave are adjusted to fit the rotational 
character of the underlying flow. Note in (3.4) that while the horizontal and 
vertical velocities are given by straightforward expressions, the complexity of 
the pressure has drastically increased compared to the irrotational case [17,18]. 
Remember that this solution is a small disturbance of the original shear flow, 
according to (3.1). For small e, we thus have an approximate solution to (2.3). 

To normalize the reference frame Stokes made a now commonly accepted pro- 
posal. In the case of irrotational flow he required that the horizontal velocity 
should have a vanishing mean over a period. Stokes' definition of the wave 
speed unfortunately cannot be directly translated to waves with vorticity (see 
[24]). In the setting of waves with vorticity we propose the requirement 

[\(x,0)dx = 0, (3.6) 
Jo 

a "Stokes' condition" at the bottom. This is consistent with deep-water waves 
(cf. [39]), and for U{y; s) it results in s = 0. As we shall see in subsection 3.1 
this indeed seems to be the natural choice of s, since this and only this choice 
recovers the well established bound y/gh for the wave speed. This is also the 
choice made in [22]. We emphazise that (3.6) is only a convention for fixing 
the reference frame; without such a reference it is however meaningless to e.g. 
discuss whether physical particle paths are closed or not. 

The corresponding approximation to the original system (2.1) is 

\{t, X, y) = -ujy + ^±g^ cos {kx - ft) cosh (ky) , 
<t^^^y) -'^W^^^(kx-ft)smh{ky), 
Pit,x,y) ^ Po + g{h - y) + ^^^^ cos (kx- ft) (3.7) 

X ^(/ -I- kuy) cosh (ky) — uj sinh (ky) 
ri{t,x) — h + acos{kx — ft) . 



Here 



, 27r ^ ^ 27rc 
k = — and / 



A ■' A 

are the wave number and the frequency, respectively. The size of the distur- 
bance is proportional to a in the whole quadruple {ri,u,v,p), so this solution 
satisfies the exact equation with an error which is O(a^) as a — > 0. Concern- 
ing the uniform validity of the approximation procedure, leading to the linear 
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system, a closer look at the asymptotic expression indicates that this solution 
is uniformly valid for 

— oo < X — ct < oo as £— >0, 

while for the vorticity we have uniform validity in the region 

SLo — 0(1) as £ — > 0. 

A rigorous confirmation of this requires a detailed analysis similar to that 
presented in [40,41], but is outside the scope of our paper. 



3.1 The dispersion relation 



The identity (3.5) can be stated in the physical variables as the dispersion 
relation 

c — s^Jgh + hu} — ^ tanh (kh) ± \/ Agk tanh (kh) + uP' tanh^ (^^)^ , (3.8) 

valid for linearized small amplitude gravity waves on a sheared current of 
constant vorticity. Note that sy^gh — hu is the surface velocity of the trivial 
solution U{y;s) stated in the physical variables. The equation (3.8) is the 
general version of the dispersion relation presented in [24, Section 3.3]. The 
authors consider waves of wavelength 27r, whence k — 1. They also require that 
u < c, and that the relative mass flux is held constant along the bifurcation 
curve for which the linearization is the first approximation. They found the 
dispersion relation 

c — = ^ ^cf tanh(/i) + ^JAgteLnh{h) + a;2tanh^(/i) 

where -Uq is the surface velocity of the trivial solution. In the more general 
case of (3.8) the problem to uniquely determine c from k, h, and uj is related 
to the fact that the requirement m < c is necessary for the theory developed 
in [24], while in our linear theory, cu and s can be chosen as to violate that 
assumption. E.g., when s — and hu < —c it is easy to see from (3.7) that for 
waves of small amplitude a << 1 the horizontal velocity u exceeds the speed of 
the wave, at least at the surface where u ^ —hu > c. The sign in front of the 
square root depends on the sign of c — s^fgh + hu). It is immediate from (3.8) 
that this expression is bounded away from 0. Positivity corresponds to the 
case dealt with in [24] , and in that case the existence of exact solutions is well 
established. Our investigation indicates that there might also be branches of 
exact solutions fulfilling the opposite requirement u > c, and as shall be seen 
below, in that case it is possible that c is negative so that there are leftgoing 
waves on a rightgoing current. In [24] it is assumed that c > 0. 
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If c — s\/~gh + huj is positive, and the vorticity is positive as well, we get a 
uniform bound for the speed of the wave. Let 

tanh hk 



h \ h 2^' ' \ h j {2-a)uj + y?f+^ V h 

meaning that 

/ tanh kh \ r— , , r— 
c < ( + s 1 ^gh <{l + s) ^gh 




If instead c — sy/gh + hu < and a; < 0, the same argument gives that 

c> - (1 + s) Jgh. 



These calculations vindicate the choice of s = 0, since in that case we recover 
the classical critical speed y/gh. 



Another comment is here in place. In [24, Section 3.3] the authors show that 
for positive vorticity, local bifurcation from shear flows requires additional 
restrictions on the relative mass flux. Again the problem is related to the re- 
quirement that u < c, and the reason can be seen directly from their dispersion 
relation stating that 

c - s^fgh + hu > (3.9) 

As a; — > —00 it forces s — > —00 to guarantee that U (y; s) < c for all y e [0, h]. 
If on the other hand a; — > 00, the inequality (3.9) admits that s — > 00. But 
for s big enough, f/(0; s) = s^/gh > c. Since we allow also u > c there is no 
corresponding restriction for positive uj. 

To summarize, we have proved 

Theorem 3.1 For a linear gravity wave on a linear current U{y;0) = —ujy 
we have 

c 7^ —huo^ 

and the dispersion relation is given by (3.8) with s = 0, where the square root 
is positive (negative) according as c + hw is positive (negative). In particular, 
if the speed and the vorticity are of the same sign, then 



\c\ < Jgh. 
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4 The phase portraits for right-going waves 



In this section we study a cross-section of the steady solution for a right-going 
wave. This corresponds to a phase-portrait analysis of the ODE-system in 
steady variables with c > 0. Since 

{x{t),y{t)) = {u{x{t),y{t),t),v{x{t),y{t),t)) 

we find that the particle paths are described by the system 



x{t) = —Loy + A cos [kx — ft) cosh (ky) 
y{t) — Asm (kx — ft) sinh {ky) 



(4.1) 



where 

^ a{f + khu) 

~ smh{kh) ^ ■ ' 

is proportional to the small amplitude parameter a. In order to study the 
exact linearised system, let us rewrite (4.1) once more via the transformation 

x{t) ^ X{t) = kxit) - ft, yit) ^ Y{t) = ky{t), (4.3) 

yielding 

'X{t) = Ak cos (X) cosh {Y)-u;Y-f 
Y{t) ^ Aksm{X)smh{Y) ^ ' ^ 

Remember that the obtained wave is a perturbation of amplitude size, and 
thus the constant A (which includes a) should always be considered very 
small in relation to u and /. Changing sign of A corresponds to the mapping 
X ^ X + TT, so we might as well consider ^4 > 0. Since we now study only 
right-going waves for which c > 0, for positive vorticity A will always be 
positive by (4.2). For large enough negative vorticity, — cj > c/h, the original 
A is however negative, meaning that the phase portrait will be translated by 
71 in the horizontal direction. This is important for the following reason: the 
presumed surface 

h + acos{X) 

attains its maximum at X — 0. Thus the crest for c + hcu > is at X — in 
our phase portraits, but at X — tt for c + hcu < 0. 



4-1 The case of positive vorticity 



Lemma 4.1 The phase portrait for the irrotational case is given by Figure 1, 
where the physically realistic wave corresponds to the area of bounded trajec- 
tories. 
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Remark 4.2 The details of this are given in [18] and a similar investigation 
is pursued in [17]. We therefore give only the main phase plane arguments for 
Figure 1. Analytic details can be found in the just mentioned papers. 



PROOF. Symmetry and periodicity of (4.4) allow for considering only Q = 
[0,7i] X [0, oo). In this region the 0-isocline, F = 0, is given by the boundary 
dn, i.e. X = 0, X = vr, and F = 0. Within holds F > 0. The oo-isocline, 
X = 0, is the graph of a smooth and convex function 

with 7(X) — >■ cxo as X y 7c/2. Here and elsewhere in this paper cosh~^ denotes 
the positive branch of the pre-image of cosh. We have 

X(X, F) > exactly when F > 7(X), 

whence X < for X G [7r/2, vr] as well as below 7(X). 

The only critical point is thus given by P = (^0, cosh^^ {f /Ak)^. Any trajec- 
tory intersecting X = below P can be followed backwards in time below 
7(X) until it reaches X = vr. For any trajectory intersecting 7(X) the same 
argument holds. Hence there exists a separatrix separating the two different 
types of trajectories, and connecting X = n with P. 
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Any trajectory intersecting X = above P can be followed forward in time 
above 7(X) and is thus unbounded. Any trajectory intersecting 7(X) can in 
the same way be followed forward in time above 7(X) and is likewise un- 
bounded. There thus exists a second separatrix, unbounded as well, going 
out from P above ^{X) and separating the trajectories intersecting X — 
from those intersecting 7(X). By mirror symmetry around X = the ciritical 
point P must be a saddle point, and the phase portrait is complete. The last 
proposition of Lemma 4.1 is the only reasonable physical interpretation. 

Theorem 4.3 For positive vorticity, uj > the properties of the phase por- 
trait are the same as for the irrotational case, cu — 0. 

PROOF. The proof is based on what we call the comparison principle, i.e. 
by comparing the phase portrait for a; > with that for a; = 0. Now changing 
CO does not affect the 0-isochnes. The change of X induced by adding the term 
coY is 

Xoj>o < Xi^=Q, (4.5) 

at any fixed point in the phase plane with Y > (where the subsripts denote 
the two different phase-portraits). Hence the velocity field is conserved wher- 
ever X < in the portrait for a; = 0, and we need only check what happens 
with the oo-isochne (which encloses all the points where X > 0). 

For any fixed X e (— 7r/2,7r/2) and a; > 0, the function 

ip(Y) ^ Ak cos XcoshY -uY - f, Y > 0, (4.6) 

is convex, satisfying ^p{0) < and ^p{Y) ^ oo as F oo, whence it has a 
exactly one zero in (0, oo). It is moreover decreasing in uj, so that if a; increases 
the solution Y of ^(Y) = increases. This means that the oo-isocline for 
a; > remains practically the same as in the irrotational case: it is a convex 
graph lying above the one for o; = 0. Just as before there is no oo-isocline for 
X e (7r/2,7r) since there ip{Y) < 0. 

Remark 4.4 We remark that according to (4.5) the wave flattens out as uj 
increases. In view of the scale X — k{x — ct) this is the same as saying that 
large positive vorticities allow only for large wavelengths. 

4.2 The case of negative vorticity 

Theorem 4.5 For negative vorticity and small amplitude a << 1 the prop- 
erties of the phase portrait are given by Figure 2. For huj > —c the crest is at 
X — 0, while for huj < —c the crest is at X — n. In particular, the steady wave 
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for hio < —c contains from bottom and up: an internal wave propagating left- 
wards, a vortex enclosed by two critical layers, and a surface wave propagating 
rightwards. 

Remark 4.6 In all essential parts this resembles the Kelvin-Stuart cat's-eye 
flow, which is a particular steady solution of the two-dimensional Euler equa- 
tions [36, Ex 2.4]- It arises when studying strong shear layers (which in our 
case means large constant negative vorticity). 
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Figure 2. The phase portrait for negative vorticity. 



In order to handle this we need to investigate the cxD-isochne for > 0. Recall 
that A = A{a) depends linearly on the amplitude (see (4.2)). 

Lemma 4.7 For negative vorticity u < 0, if a> is small enough so that 

^=i„h-'(^)-J^77^-^ (4.7) 

a \aj V Va/ a 

is positive for a = Ak, then the 00-isocline of (4.4) for X G [0, vr] consists of 
two disjoint parts: 

(1) the graph of an increasing function Yi{X) defined for X G [0,7r], and 

(2) the graph of a decreasing function Y2{X) defined in (7r/2,7r]. 

We have Yi{X) < Y2{X) — > 00 as X \ 7r/2, and for any 5 > there exist 
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Y* > and a possible smaller a such that the slope satisfies 




m 



R= [0,7r] X [Y*,Y* + S]. 



PROOF. Just as before (p{Y) as in (4.6) is convex for X e (— 7r/2, 7r/2) with 
</?(0) < 0. However, as 



along the oo-isochne ip{Y) — 0. According to the ImpUcit Function Theorem 
[42, Theorem 1. 1.1] the curve can be continued across this point into X G 
(7r/2, 7r]. There cosX < 0, and consequently (p{Y) is now concave with ip{0) < 
0, <yj(0) > 0, and ^{Y) — > — oo as y — > oo. The function (p{Y) attains its 
global maximum when 



Thus the equation (p(Y) = has no, one, or two solutions according as (4.7) 
is negative, vanishing, or positive, for a = — Ak cos X. 

It is easy to see that if a > is small enough this expression is positive, 
while it becomes negative for large a. In view of that a vanishes as X \ 
7t/2, we see that at X = 7r/2 a new branch of the oo-isocline appears from 
Y — +00. Keeping in mind that (p{Y) is concave for X e (7r/2, tt), where cos X 
is decreasing, it follows that the upper branch of ip{Y;X) = is decreasing 
as a parametrization Y{X) while the lower branch is increasing in the same 
manner. 

Depending on the relation between A, cu, and /, it may be that the two 
branches both reach X = tt separately, that they unite exactly there, or that 
they unite for some X < tt, where they cease to exist. However, if A is small 
enough in relation to |a;| and /, (4.7) guarantees that both branches of the 
oo-isochne exist as individual curves throughout X e (7r/2, vr]. 

For the final assertion, remember that the slope is given by 



Fix Y* with -ujY* > {1 + f + 6). Since A ^ as a ^ there exists oq such 
that for any a < ao the inequality Ak cosh.{Y* + 5) < S/n holds. In view of 
that sinh ^ < cosh ^ this proves the lemma. 

We are now ready to give the proof of Theorem 4.5. 



X n/2 we now have Y 





dY .4/,- sin X sinh y 



(4.8) 



dX Ak cos X coshY -LoY - f 
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PROOF. By periodicity and horizontal mirror symmetry it is enough to 
consider fl = [0, tt] x [0, oo) (remember that F = is the bed). The first critial 
point is 

Po = {0, Y) where Ak cosh Y - uuY - f = 0, 
while the second and third critical points are 

Pi = (7r,Yi) and P2 = (7r,Y2), 

where Yi and Y2 are, in order of appearance, the smallest and largest solutions 
of 

Ak coshY + u;Y + f ^0. 

In n holds y > 0, while the sign of X is negative below Yi{X) and above 
Y2{X), respectively positive elsewhere in Q. This follows from Lemma 4.7, 
and can be confirmed by considering dyX for a fixed X. 

The system (4.4) admits a Hamitonian, 

H{X, Y) = Ak cos X sinh Y-^ uY^ - fY, (4.9) 

with 

X = dyE, 

Y = -dxB, 

for which the trajectories of (4.4) are level curves. To determine the nature of 
the critical points we study the Hessian, 

(cos X sinh Y sin X cosh Y 
sin X cosh Y — cos X sinh Y + — — 
Ak' 

At Pq) where X = 0, we immediately get that there exists one positive and 
one negative eigenvalue, whence the Morse lemma [43] guarantees that Pq is 
a saddle point. Insertion of P2 = (7r,y2(7r)) yields 

/sinhr2(7r) \ 
D'HiP2) = Ak\ . n ,w X ^ • (4.10) 

Now remember that P2 is the point where the function v^(l^), for X = tt, at- 
tains its second zero. This happens when the derivative '^'(Y) — — AA;sinhy — 
a; < 0, and consequently also P2 is a saddle point. 

That Pi is a center can be seen in the following way: D^H{Pi) differs from 
(4.10) only in that Y2{n) is substituted for Yi^n). Since Pi is the point of the 
first zero of (p{Y) for X — n^it follows that there (p'{Y) = —Ak sinh F — a; > 0, 
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whence the Hessian is a diagonal positive definite matrix. According to the 
Morse lemma there exists a chart {x, : — > such that 

H{X, Y) = H{P^) + x\X, Y) + y\X, Y) 

in a neighbourhood of Pi. (An alternative and efficient way is a phase por- 
trait argument using the symmetry, which ensures that any trajectory that 
intersects X = n twice is closed.) 

Pq being a saddle point, there is a separatrix 7f - i.e. a trajectory separating 
two qualitatively different trajectory behaviour^ - which can be followed 
backwards in time from Pq below Yi{X), and a separatrix 7^ which can be 
followed forward in time from Pq above Yi{X). By the direction of the velocity 
field, 7f connects Pq with X = n below Pi. 

In the same manner there are separatrices 7^ and 74*" leaving the saddle point 
P2, and since 74*" lies above Y2{X) it is unbounded and encloses a family of 
unbounded trajectories starting from X = tt above P2. The separatrix 7^ can 
be followed backward from P2 below Y2{X). 

Now, according to the last part of Lemma 4.7, there is a family of trajectories, 
{P}, starting from X = above Pq and reaching X = vr in finite time in 
between Pi and P2. Following 7^ we therefore must intersect X = above 
{P}. In general, starting from Y2{X), following any trajectory backwards, we 
find that it must intersect X = above {P}. Since X\y2{x) = the same 
trajectory starting from Y2{X) is contained above Y2{X) and is unbounded. 
The family {P} also guarantees that 72" connects Pq with X = tt above Pi 
but below P2. The phase portrait is thus complete. 

Remark 4.8 In case (4.7) does not hold to be positive one might still pur- 
sue the analysis, finding a fluid region where the situation is the same as for 
positive and vanishing vorticity. However, above the surface the situation is 
radically different and, for fixed amplitude and other parameters, indicates a 
transition between negative and positive vorticity. In particular, there exists an 
a; < for which a bifurcation takes place: a second critical point appears, and 
as uj decreases it immediately gives birth to a third critical point. Since our 
linearized model presupposes that the amplitude is small, we do not investigate 
this transition further. Some of the main features are given, without proof, in 
Figure 3 below. 



The use of the word separatrix is somewhat ambigous. In our case, however, the 
geometrical definition corresponds to the analytic notion of the stable and unstable 
manifolds which are defined by H(X,Y) = H(Pi), i = 0,2, and whose existence 
follow from by the Implicit Function Theorem [42, Theorem 1. 1.1]. 
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Figure 3. The bifurcation of the phase portrait. 



5 The physical particle paths 

In this section we shall investigate how from the behaviour of the trajectories 
{X{t),Y{t)) we might infer the motion of the physical particles {x{t),y{t)). 
Remember that the relation between those two pairs are given by (4.3). 

Before moving on we recall from [17,18] that in the case of irrotational linear 
waves, the particles were found to move in almost closed orbits, with a slight 
but positive forward drift. This is in line with the classical Stokes drift [12] 
according to which there is a forward mass drift. It is also consistent with 
results for the exact equations [20,21], showing that within regular Stokes 
waves - which are irrotational - no orbits are closed. 

We emphasize that the appearance of vorticity radically changes the picture. 
In particular, it will be shown that both when the vorticity is negative and 
satisfies a + hu < 0, and in the case of large positive vorticity, there does not 
exist a single pattern for all the fluid particles. Rather, different layers of the 
fluid behave in qualitatively different ways, with some layers moving constantly 
in one direction. For the same reason, it is hard to state any transparent results 
other than Theorem 5.3 stating that for small amplitude waves on a current of 
large positive vorticity there are indeed closed orbits, and Theorem 5.5 which 
asserts that for negative vorticity all fluid particles display a forward drift. 
Figure 4 however shows the main features for negative vorticity, as well as a 
possible situation when the vorticity is positive and very large. 

Lemma 5.1 Particles near the flat bed y = have a forward (rightward) 
drift. 



PROOF. Consider the time r that it takes for a particle {X(t),0) with 
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X{0) — TT to reach X{t) — — tt. We have 

dt{X) ,^ dX 



f -Ak cosX 



Tr/2 \ f — Ak cos X f + Ak COS X 

/2 f^{Ak COS xy " f 



r/^ dX 2n 



We assert that this holds also near the bed: for any fixed X we may differen- 
tiate X{X, Y) with respect to Y, obtaining 

dyX — Ak cos X sinh Y — u. 

By continuity, there exists S{e) > such that \X(X, Y) — X{X, 0)\ < e when- 
ever < Y < 6, uniformly for X G M. If we thus consider r for a trajectory 
intersecting X = at level Y e (0, 6) we may choose e arbitrarily small so to 
obtain r > 27r//. 

Now a closed physical trajectory implies y{T) = y{0) for some T > so that 
Y{T) = Y{0) in view of (4.3). It follows from the phase portraits that for 
trajectories close enough to the bed, this forces 

X{T) - X{0) = -2nn, meaning = x{T) - x{0) ^ fT - 2nn, (5.1) 

for some n G N. By periodicity, T/n = r is the time it takes the trajectory 
X{t) to pass from X — tt to X — — tt (which any trajectory near the bottom 
does). Prom (5.1) we infer that T — 27r/f. 

Remark 5.2 We also see from this reasoning that if r > 27:/ f , then the 
particle will be to right of its original position, and contrariwise. 



5.1 The case of positive vorticity 



Theorem 5.3 If the vorticity oo > is large enough, and the amplitude small 
enough, then there particles moving constantly forward as well as particles 
moving constantly backward. In particular there are closed orbits. 



PROOF. By Lemma 5.1 the particles nearby the flat bed y = display a 
slight forward drift. In principle, they behave as in the irrotational case (see 
[17,18]). 
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By continuity X(X, Y) can be made arbitrarily small, uniformly for all {X, Y), 
close enough to the critical point Pq. Thus for the trajectories near Pq the time 
r as in Lemma 5.1 can be made arbitraily large, and hence there is a forward 
drift i; > for the corresponding physical particles. 

In between those two layers something different might happen. Fix Y*, 5 > 0, 
and choose < a <C 1 small enough such that 

Akcosh(Y* + 5) < 5. 

Then choose a; ^ 1 such that ujY* — 6 > ti, and such that the solution Yq of 
Ak coshYo — cuYq — / = satisfies Yq > Y* + 6 (cf. (4.6) and the paragraph 
following it). Then the slope given by (4.8) satisfies 

\dY/dX\ < S/n, 

so that the trajectory for which 1^(0) = Y* + 6 remains in [Y*, Y* + 6] where 
X < — TT — /. Hence 

x(t) < 0, 

for all t for the physical particle and there is a constant backward drift. Since 
the physical surface is given hy Y = k{h + acosX) we can adjust our choices 
so that Y* + 5 < k{h — a) guarantees that the orbits we consider are indeed 
within the fluid domain. 

Using continuity once more we find that for large positive vorticity and small 
amplitude, there do exist closed physical orbits. 

Remark 5.4 Since the physical crest might lie below the critical point Pq we 
can only be sure that there is at least one (infinitessimally thin) layer of closed 
orbits. However, this appears to happen only for small a and large uj, and even 
so it does not have to affect more than one single trajectory in the {X, Y) - 
plane. The particle paths are depicted to the left in Figure 4- 

5.2 The case of negative vorticity 

In the case of irrotational linear waves, it was found in [17,18] that all fluid 
particles display a forward drift. This is confirmed for negative vorticity, and 
the proof of Theorem 5.5 also shows the situation in the three different layers 
of the fiuid. A schematic picture of this can be found to the right in Figure 4. 
Recall that for c + /ic<j > all of the fluid domain lies beneath the lowest 
separatrix, so that the situation is the same as for zero vorticity. For c+huj < 
the fluid domain streches above the vortex so that the picture is quite different 
from irrotational waves. 
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Theorem 5.5 For negative vorticity, all the fluid particles have a forward 
drift. 



PROOF. For reference, consider Figure 2. We treat separately 

i) the interior wave (beneath the lowest separatrix), 

ii) the vortex (between the first and second separatricies from bottom and up), 
and 

iii) the surface wave (between the second and third separatricies from bottom 



The case i). Any trajectory {X{t), Y{t)) in this region passes X = kn, k G Z. 
We may thus consider r = ^ dX as in Lemma 5.1, r being the time it 
takes for the particle to travel from X{0) — n to X{t) = — tt. Again, for any 
fixed X, 

dyX = Ak cos X sinh Y - uj > -uj - Ak sinh Y = (p'{Y) 

with the notation of (4.6). Since </?' has its only zero at the level of Pi, it 
follows that dyX > in the interior wave. Since X < at y = 0, we deduce 



so that - according to Remark 5.2 - the physical particle path describes a 
forward motion. 

The case ii). Any trajectory {X{t),Y{t)) within the vortex is bounded and 
passes -'^(O) = tt, whence 



and up). 



that 




x{t) = x{0) + 



ft + B{t) 



where \B{t)\ < tt. 



k 



In particular, at Pi we have X = 0, so that the physical particle moves straight 
forward according to x{t) — {n + ft)/k, for alH > 0. 



The case iii). We need only observe that whenever X is positive so is i; = 
{X + f)/k, whence all trajectories above the 0-isochne connecting Pq and Pi 
correspond to fiuid particles moving constantly forward. 
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Figure 4. To the left the physical particle paths for small-amplitude waves with 
large positive vorticity is depicted. The black arrows describe what happens to 
some typical particles in time, while the axis marked 5 correspond to the sepa- 
ratrix of Figure 1, and the C is the critical point in the same figure. Note that 
depending on the amplitude and the vorticity, the surface of the physical wave need 
not correspond to the uppermost arrow. Theorem 5.3 however guarantees that for 
large enough vorticity and small enough amplitude the surface lies strictly above 
the closed particle path (the first circle from bottom and up), so that there are par- 
ticles with a mean backward drift as well as the opposite. To the right we see the 
particle paths when the vorticity is negative such that c -|- huj < 0. This corresponds 
to Figure 2 with the same notation as above. Near the bottom we have a mean 
forward drift with nearly closed ellipses, but within the vortex of Figure 2 we see 
a drastic change of behaviour with a constant forward drift. This is retained even 
above the separatrix separating the vortex from the surface wave. 




6 Summary and discussion 



We have deduced and investigated the closed solutions of linear gravity water 
waves on a linearly sheared current (constant vorticity). Such linear waves 
satisfy the exact governing equations with an error of magnitude a^, where a 
is the amplitude of the wave. The main purpose has been to understand how 
the presence of vorticity influences the particle paths. While in the irrotational 
case all the particles describe nearly closed ellipses with a slight forward drift, 
we have found that vorticity might change the picture. For positive vorticity 
the situation is very much the same as in the irrotational case, but for large 
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enough vorticity and small cnoTigh waves, there are closed orbits within the 
fluid domain. For negative vorticity exceeding the wave speed sufliciently much 
all the particles describe a forward drift, but the nearly closed ellipses can be 
found only in an interior wave near the flat bed. 

It seems that all waves of constant vorticity are qualitatively though not quan- 
titatively the same unless we accept the speed of individual particle to exceed 
the speed of the wave. Then appears waves with interior vortices. So far there 
is no corresponding exact theory of such rotational waves, since all work has 
focused on regular waves not near breaking and without stagnation points. 

When discussing particle paths it is important to remember that the question 
of closed orbits is valid in relation to some reference speed. For irrotational 
waves Stokes required that the average horizontal velocity should vanish. For 
waves with vorticity we propose that the same requirement at the bottom is 
the most sensible counterpart of Stokes' deflnition. This is supported by the 
fact that only for that choice we recover the classical critical wave speed \fgh. 

While interesting in its own right, the investigation pursued here might have 
further implications for the numerous and well-known model equations for 
water waves, e.g. the Kortweg-deVries, Camassa-Holm, and Benjamin-Bona- 
Mahony equations. They all describe the surface - or nearly so - of the wave. 
Though reasonable for irrotational waves, flndings on uniqueness for rotational 
waves indicate the same as our investigation: beneath two identical surfaces 
there might be considerable different fluid motions (see Figures 1 and 2). Apart 
from the trivial case of a flat surface there are so far no known exact examples 
of this possible phenomenon, but if true it might motivate a new understanding 
of in what sense the established model equations model the fluid behaviour. 
Indeed vorticity, even when constant, is a major determining factor of the fluid 
motion, and it should as such be considered highly important in the study of 
water waves. 
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